Collective oscillations of dipolar Bose-Einstein condensates and 
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We propose a scheme for the measurement of the s-wave scattering length a of an atom or molecule 
with significant dipole-dipole interaction with an accuracy at the percent level. The frequencies of the 
collective oscillations of a Bose-Einstein condensate are shifted by the magnetic dipole interaction. 
The shift is polarization dependent and proportional to the ratio £dd of dipolar and s-wave coupling 
constants. Measuring the differences in the frequencies for different polarization we can extract 
the value of £dd and thus measure a. We calculate the frequency shifts for a large variety of non- 
axisymmetric harmonic traps in the Thomas-Fermi limit and find optimal trapping geometries to 
maximize the shifts. 
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Recent research on polar molecules Q , Rydberg atoms 
3], laser-induced dipoles dipolar spinor condensate 
$ and Bose-Einstein condensates (BECs) of polarized 
chromium [^, Q is attracting a growing attention towards 
the exciting new physics of dipolar gases. Dipolar gases 
have intricate stability diagrams and ground state prop- 
erties 0, Ijj fi cill , including new quantum phases in op- 
tical lattices[ri]| and roton dispersion relations 

In this paper, we present exact semi-analytical results 
for the frequencies of the collective condensate oscilla- 
tions in non-axisymmetric traps. We use recent analyt- 
ical results [fj, |!J obtained in the experimentally rel- 
evant Thomas-Fermi (TF) regime of a large number of 
atoms or molecules N . Previous studies of collective ex- 
citations were confined to the case of axisymmetric traps 

SEEi. 

As an application of our results, we discuss a new 
method for the accurate determination of the s-wave scat- 
tering length a at the percent level. This possibility is 
entirely due to the dipolar interaction and its high degree 
of manipulability. Measurements of a at 1% accuracy 
level near a Feshbach resonance may be used to detect 
ultra-small variations of the fundamental constants |l5| . 

The density n of a dipolar BEC and its superfluid ve- 
locity v obey the continuity and Euler equations given 
by ' 
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where Vho = ( TO /2) (u 2 x 2 + uo 2 y 2 + w 2 z 2 ) is the trapping 
potential and V m { is the mean field potential given by 
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where m is the mass of the dipolar particle. Vdd is the 
anisotropic and long-range dipolar interaction given by 
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where Cdd = MoMm f° r magnetic particles of moment \i m 
and Cdd = d 2 /eo for particles with an electric dipole 
moment d e . The properties of dipolar BEC depend ad- 
ditionally on the dipole orientation and on the dimcn- 
sionless ratio between the dipolar coupling Cdd and a 
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An exact class of solutions of the TF hydrodynamic 
equations and @ has been obtained in 0, 0, 0] . 
The existence of this class of solutions for harmonic traps 
is due to the harmonic nature of the external and self- 
consistent potentials and the Bernoulli term. The density 
has the "usual" form of an inverted parabola given by 
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The velocity field v{r, t) = (1/2) V [a x x 2 + a y y 2 + a z z 2 ] 
depend on the condensate radii Rj with j = x,y,z 
through the relation atj = dt \og(Rj). The Rj are ob- 
tained by solving the Newton equation 

Rj = —dju(R x ,R y ,R z ) 



(7) 
The 



where dj denotes the partial derivative d/du, 
reduced potential landscape u is the sum of Uho = 
(1/2) (w 2 i? 2 + oj 2 R 2 + uj 2 R 2 z ) from the harmonic trap, 
u s = (lhNh 2 a/m 2 ^ / R x R y R z from the s-wave pseudopo- 
tential and Udd = — £dd K2)u s from the dipole-dipole 
interaction. For dipoles aligned with the x axes the con- 
densate aspect ratios ki,K2 are given by k\ = R y /R x 
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and K2 = R z /R x . The symmetric function / with values 
in the interval (—2, 1) is given by 

- H3«,.„m^ (8) 



with sin ip = (1 — k 



i 1 / 2 and sin 2 a = (1 — k 
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F(ip\a) and E(y\a) are the incomplete elliptic integrals 
of the first and second kinds (see |6j for more details). 

The equations for the small oscillations around equi- 
librium are obtained expanding the Newton equation (JZJ) 
for the condensate radii Rj around the equilibrium values 
R £ a for small displacements SRj 



Rj - R e j 



SRi 
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where 



u" m = dl m u\ Ri 



(9) 



(10) 



is the matrix of the second derivatives of u with respect 
to the condensate radii Ri and R m evaluated on the equi- 
librium values The frequencies 0^ of the small oscil- 
lations around equilibrium are simply given by the square 
roots of the eigenvalues of u" m . The index k of the normal 
modes are chosen so that < f2 2 < 

Fig. 1 shows an example of the dipole-induced shifts 
on the frequency of the collective modes as a function of 
£dd between and 1 for a possible experimental relevant 
anisotropic trap. Cdd can be varied by tu ning the scatter- 
ing length via a Feshbach resonance (see [lj| for the case 
of chromium). Fig. 1 shows how for sufficiently low val- 
ues of Edd the f2fc depend almost linearly on edd- In the 
following, we shall employ this fact to obtain simplified 
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Figure 1: Solid, dashed and dotted lines represent in (a) 
and (c) the three collective frequencies Q3, f^2 and f2i (in 
units of uiz) and in (b) and (d) the corresponding relative 
variations AO/fi as function of £dd in the TF interval (0, 1). 
The magnetic dipoles are polarized along the i-axes in (a) and 
(b) and along the y-axes in (c) and (d). The trap frequencies 
are u x = 27r942 Hz, uj v = 27r712 Hz and uj z = 27iT28 Hz. 



expressions which allow for a quick and simple estimation 
of the dipole-induced shifts. 

In the absence of dipole-dipole interaction the unper- 
turbed matrix w" m is given by the simple expression 

u'{ m \ = 2ujfSi m +u!iuj m (11) 

where |o denotes the evaluation at edd = 0. The eigen- 
values and eigenvectors of (|11|) can be calculated analyt- 
ically but their expressions are too long to be reported 
here. To first order in edd the variation of the matrix u[' m 
is given by 



Hm 
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where di m Udd is the dipolar contribution to v!( evalu- 
ated on the unperturbed equilibrium configuration i?^|o 
and 



{Rf dijiu)\ = 



- 2ojf6. 



2ujj Sij 



uji uu l 5ji - lui ujj 5 a - uji cuj (13) 



is part of the contribution 9/ nm u|o AR^ due to the dipo- 
lar displacement of the equilibrium Ai?^|o = R e n — -R^Jo 
given by [(| 
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where / and its partial derivatives d Kl /, d K2 f are evalu- 
ated in Ki = uj x /ujy and k 2 = Lu x /uj y for dipolc polarized 
in the x-direction. Using standard first order perturba- 
tion theory in edd we obtain for the relative frequency 
shifts 



2(n h 



(17) 



where are the corresponding normalized eigenvectors 
of the unperturbed matrix (|11[> . 

Fig. 2 shows the results of the linear approximation in 
edd- These results allow a quick estimation of the dipo- 
lar frequency shift of the two lower collective modes for 
different geometries. The upper mode (monopole-like) 
frequency shifts are not shown in Fig. 2 as they repre- 
sent a weaker effect of about a factor two. As most of 
the traps have at least two frequencies that do not differ 
significantly, we consider the cylindrical symmetry as a 
starting point (the fifth row in Fig. 2) and we modify the 
two originally equal trap frequencies by factors ranging 
between 4 and 0.35. We observe that non-axisymmetric 
traps maximize the frequency difference for different po- 
larization normalized by the unperturbed frequency. In 
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Figure 2: Dipolar relative frequency shifts of the collective modes obtained after linear expansion in £dd (for convenience 
£dd = 0.15). These shifts depends both on the trap aspect ratios \ v = u x /ujy, \ z = u) x /uj z and dipole polarization. The first 
two columns show the relative frequency shifts AQ/Q as a function of the logarithm of the trap aspect ratio X y . The column 
on the right shows the Log-Log plot of the unperturbed frequencies as a function X y for reference. Dashed lines and solid lines 
correspond to the intermediate and lower frequency modes, respectively. Different rows of figures correspond to different ratios 
of Xz/Xy and polarization. Figures on the left (center) have polarization B\\x (B\\y). From the top \ Z /X y = 4, 2.83, 2, 1.41, 
1.00, 0.71, 0.50, 0.35. Note in the case A z = X y , which corresponds to cylindrical symmetry, the dashed and solid lines have a 
discontinuity in X y = 1 (that corresponds to the spherical symmetry), due to a level crossing of the two quadrupolar modes. 
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particular, the relative frequency difference for the in- 
termediate mode reaches a value as large as 0.057 for 
2. 



A 2 > lAIXy and X y 




210 220 230 240 250 260 
t [ms] 



Figure 3: A modulation of 5% of the harmonic trapping 
potential at frequency of about Q2 = 2-7T950 Hz is applied for 
5 oscillation periods. The R y radius (arbitrary units) of a 
chromium BEC (edd ~ 0.15 ± 0.2) exhibits large oscillations 
in time of about 50%. Solid line and dashed line correspond 
to the x and y polarization. The trap frequencies are assumed 



2tt1000 Hz, 



= 2tt600 Hz and lu z = 2tt200 Hz. The 



modulation frequency is chosen so that the second mode is 
excited. The phase shift is quite visible already from the 
beginning and after 230ms (~ 34 — 36 periods for the two 
polarizations) is about 2 x 2tt. 

Finally, we discuss how from the measurement of the 
oscillation frequencies of a condensate for different polar- 
izations we might determine a to the percent level. As 
an example, we consider a chromium BEC confined in 
a trap with quite distinct frequencies lo x = 27rl000 Hz, 
u) v = 27r600 Hz and uj z = 27r200 Hz. In this case, it is pos- 
sible to excite each collective mode just modulating the 
laser intensities at approximatively its mode frequency 
for a few periods of oscillations. Fig. 3 illustrates the 
process of excitations of the second mode which has a 
frequency of about CI2 = 27r950 Hz. The mode oscilla- 
tion in mainly concentrated along the y-axes (the axes of 
the middle trap frequency). Fig. 3 shows the time evo- 
lutions of the condensate radius R y for the two different 
polarization by applying a modulation of 5% of the laser 
intensities for 5 oscillation periods T. We assume that we 
can see 30-40 clear oscillations. When the approximate 
values of each mode frequency can be inferred from other 
data we do not need to sample the full time interval and 
the measurements can be concentrated only at the begin- 
ning (t ~ 35 ms) and around t ~ 230ms (about 35 free 
oscillations). We estimate that we can measure each fre- 
quencies at least with an accuracy of 10 -3 . In the case of 
chromium (edd ~ 0.15), this translates into an accuracy 
of about 2 x 10~ 2 for the scattering length, which would 
improve by one order of magnitude the 20% accuracy of 
the today's most recent measurements [l7j |. 
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